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Abstract 

The distribution of energy level separations for lattices of sizes up to 
28x28x28 sites is numerically calculated for the Anderson model. The re- 
sults show one-parameter scaling. The size-independent universality of the 
critical level spacing distribution allows to detect with high precision the crit- 
ical disorder W c = 16.35. The scaling properties yield the critical exponent, 

v = 1.45 ± 0.08, and the disorder dependence of the correlation length. 
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The statistical properties of the energy spectra of disordered quantum systems in the 
metallic regime can be described by the random-matrix theory initiated by Wigner and 
Dyson This was shown by analytical approaches such as the super symmetric nonlinear 

cr-model formalism Q and the impurity diagram technique ||. A characteristic property in 
this regime is the repulsion between adjacent energy levels which is caused by a pronounced 
overlap of delocalized one-electron states. 

By increasing the fluctuations of the random potential in three-dimensional (3D) systems 
a delocalization-localization transition is induced, the Anderson metal-insulator transition 
(MIT) ||. In the insulating phase, when the magnitude of the fluctuations of the potential is 
much larger than the critical value, the states of the system are localized. The corresponding 
energy levels are not correlated due to the vanishing of their spatial overlap. Therefore, 
spectral fluctuations are governed by the statistics of independent random variables. A 
central question in discussion of the level statistics in disordered systems is how the spectral 
fluctuations change their behavior when extended states are transformed into localized ones 
when the disorder is increased. 

According to the one-parameter scaling theory ]7| the conductance of the system shows 
critical behavior at the MIT. Since the conductance is related to the spectral two-level 
correlation function one can expect that the level statistics also experiences the crossover 
between the metallic and the insulating regimes |§. Correlations of the energy levels are 
contained in the distribution of the distances between successive energy levels, P{s). In the 
metallic region P(s) can be well described by the Wigner surmise 

7T f 7T \ 

p w(s) = -s exp [~^s 2 j , (1) 

proposed primarily for the Gaussian orthogonal ensemble of random matrices whose elements 
fluctuate normally around zero with the equal variances Here the distances between 

neighboring energy levels s are measured in units of the mean level spacing A. In the 
localized regime the spacings are distributed according to a Poisson law 

P P {s) = exp(-s). (2) 
2 



The crossover from the metallic to the insulating regime is accompanied by a change in the 
level spacing distribution from correlated to completely uncorrelated behavior. One can, 
therefore, suspect that the properties of P(s) for finite systems, when investigated as a 
function of the size of the system, would allow for the detection of the MIT and even the 
critical properties of the transition. As a consequence, the crossover of P(s) between the 
Wigner ([[]) and the Poisson (§) statistics was extensively studied both analytically |9|Jl0t| and 



numerically . By using of the finite-size scaling argument for P(s) a new criterion 



for the location of the metal-insulator transition was proposed by Shklovskii et al ]TTJ . This 
argument was cross-checked by numerically determining level statistics, and investigating 
the scaling properties for large s. One advantage of this method is that one needs only the 
energy spectrum without calculating eigenfunctions or the conductivity. 

In this paper we focus again on the critical behavior of the level statistics near the 
disorder-induced MIT. We study the finite-size scaling properties P(s) in the whole range of 
spacings. For the calculation of the energy levels, numerical diagonalization of the Anderson 
model was performed by using a modified Lanczos algorithm for systems of sizes up to 
28x28x28 lattice sites. We determine very precisely the critical value of the disorder W c and 
the correlation length exponent v. By considering the integrated distribution of spacings 
in various regions of P(s), we find W c = 16.35 and v = 1.45 ± 0.08 in good agreement 



with similar earlier calculations [I3J and previous findings based on the transfer-matrix 



(TM) method for the modulus of the Green's function |15|JT6[] . We provide an independent 
determination of the correlation length as a function of the disorder, £,{W). Our results 
demonstrate unambiguously and strikingly, that the integrated level spacing distribution 
plays not only the role of a scaling variable at the Anderson transition, but can also be used 
to define qualitatively and accurately the critical properties. 

The Anderson model for a disordered system is described by the Hamiltonian H = 
J2n £ n c n c n + ^n^m( c n c m + c n c rn)i where and c n are the creation and annihilation operators 
of an electron at a site n in a lattice; m denotes the number Z of the nearest neighbors of 
the site n (Z = 6 in a 3D simple cubic lattice). The site energies e n measured in units 



of the overlap integral between adjacent sites correspond to the random potential. They 
are uniformly distributed in the interval [-W/2, W/2], where W specifies the degree of the 
disorder of the system. It was found earlier by the TM-method (15| that the critical disorder 
of the MIT which occurs in the middle of the band E = corresponds to W c ~ 16.5. Thus, 
for a finite system when W decreases from W c , the level statistics should approach ([I]), 
and when W 3> W c it should obey @ J8fl. In order to find the electron spectrum in the 
critical region for W close to W c we diagonalized numerically the real symmetric matrices 
of the Anderson Hamiltonian for simple cubic lattices of the size L 3 with periodic boundary 
conditions. One has to consider system sizes much larger than the lattice constant. This is 
necessary in order to get rid of the microscopic peculiarities of the discrete lattice model, 
and to be able to compare numerical with analytical results obtained in the thermodynamic 
limit. Despite H is represented by a very sparse matrix the computation of the energy levels 
of large systems is highly nontrivial due to the huge matrix size which grows as L 6 . 

We have developed a considerably improved Lanczos algorithm for the eigenvalue problem 
of very sparse and big matrices. One of the improvements is based on the hierarchical 
structure of H which consists of three levels of complexity depending on the dimensionality 
of the lattice. Another feature is that all non-zero off-diagonal matrix elements are equal 
to unity. This enabled us to decrease substantially the number of multiplications during 



the recursion [fL7|j . In addition, the absence of degeneracy of the eigenvalues leads to a 
faster convergence of the eigenvectors. We exploited the effect of level repulsion in the 
range of s < 0.1 for diminishing the number of steps of the Sturm sequence in the bisection 
method that is used for searching the eigenvalues of the tridiagonal matrix produced after 
the recursion. In this case one can govern the precision conditions for resolving two very close 
'real' eigenvalues provided that the absolute accuracy for recognizing 'spurious' eigenvalues 
is kept as high as of the order of 1CT 12 . All of these optimizations enabled us to create a 
computer program which consumes much less time than a conventional implementation of 
the Lanczos method |T7|] . As a result, we succeeded in diagonalizing matrices of sizes up to 



21952 x 21952. We applied the algorithm to cubes of linear sizes L = 6, 8, 12, 16, 20, 24, 28 



at various degrees of the disorder close to W c . 

One consequence of the fact that the distribution of e n is box-like is that the correlation 



length £ depends only very weakly on the energy near the band center |16j| . This follows 
from the phase diagram of localization which shows that W C (E) is almost constant when 
the mobility edge varies from E = —6 to 6. This allows us to use a wide central part of the 
band containing half of all of the eigenvalues of H, without mixing localized and delocalized 
states. One can check that the level statistics P{s) for narrower energy intervals do not 
change considerably. The calculations were carried out for ensembles of different random 
configurations. For a given pair of {L, W} the number of samples was chosen such that a 
total of approximately 10 5 eigenvalues were taken into account. Since the energy interval 
is so wide, the density of states p = (AL 3 ) -1 depends on energy, though only slightly, one 
has to unfold the spectrum. The unfolding was performed by using parabolic spline fits for 
p{E). 

As an example, Fig. [I] displays the histograms of P{s) for L = 6 and 28 at various disorder 
W near the MIT. By increasing W the spacing distributions for both L vary continuously 
from Pw{s) to -Pp(s), particularly in the range of s ^ 2. Results for other, intermediate, 
sizes confirm this smooth crossover of P{s). However, the behavior of the crossover depends 
substantially on the size of the cube. One can observe in Fig. [I] that P{s) changes from ([!]) 
to (0) faster for L = 28 than for L = 6. As the size of the system increases the function 
P{s) approaches the above two limiting regimes within a narrower interval of W around 
W c . The size- dependence of P(s) is observed on both sides of the MIT. Close to W = 16.5, 
the spacing distribution has almost the same form for all considered L, from 6 up to 28. 
The independence of P{s) on L at W = W c is consistent with the prediction of a novel 
universal level statistics at L —>■ oo that exists exactly at the transition [|TTJ. One can expect 
that in the thermodynamic limit there are three asymptotic shapes, namely, the Wigner 
surmise (|1|) for W < W c ] the Poisson law (0) for W > W c , and the critical distribution P c (s) 
when W = W C . 

Since P{s) scales differently with L for various regions of s, it is more convenient to 



consider the integrated probability J(s) = P(s') ds', which equals a fraction of those 
spacings which are larger than s. Because all spacings are positive, 1(0) = 1. In addition, 
J °° I(s) ds = 1 due to normalization to the total number of spacings in a given energy 
interval. Eqs. ([I]) and (|2|) give I\v(s) = exp(— 7rs 2 /4) and Ip(s) = exp(— s) for the two 
limiting regimes. For numerical purposes we introduce J(s) = (I(s)— Iw(s))/(Ip(s) — Iw(s)). 
It describes the relative deviation of I(s) from Iw( s )- 111 order to extract the dependence 
only on W and L and to minimize numerical errors we define a certain value of the spacing s Q 
for J(s). The choice of the average spacing, s Q = 1, is not favorable, because J(A) exhibits 
a very small size effect. One can see in Fig. [I] that Pw{s) and Pp(s) cross at s\ = 0.473 
and s 2 = 2.002. Therefore we studied a(W, L) = J(si) representative for the small-s part of 
P(s), and j(W, L) = J(s 2 ) representative for its tail, where Iw( s ) an d Ip(s) are significantly 
different. The former reflects the "strength" of the repulsion of two consecutive levels when 
the separation between them is less than s\, while the latter is the exponentially small 
probability to find a gap in the spectrum which is wider than two mean spacings A. For 
L — ► oo both, a and 7, are equal to zero for W < W c and equal to unity for W > W c . 

Fig. shows the disorder dependences of a and 7. The curves for different sizes cross 
at one point, W c , which corresponds to the disorder-induced MIT. Near the critical point, 
as long as \W — W c \ ^ 2, a(W) and j(W) depend stronger on W for larger cubes. A 
similar behavior of 7 was observed earlier 0] for much smaller cubes than ours. A different 



variable, J S2 ° I(s) ds + s 2 — 1, was studied in f["4jl , which corresponds to using P(s) in the 



range of large spacings, as 7 does. A comparison of the two plots in Fig. ^| shows that 
near the MIT the scaling variables a and 7 behave almost identically. The critical disorder 
and the scaling properties of the integrated level spacing distribution do not depend on the 
choice of s Q . We have checked this s G -independence by repeating the calculations for s o =0.1 
and 4. However, the results for s± yields the highest accuracy. Using the data for a one can 
determine very accurately the critical value of the disorder, W c = 16.35 ± 0.15. 

In the following we discuss the critical behavior of a. It characterizes the slope of the 
linearly growing part of P{s). Assuming the validity of the one-parameter scaling a(W, L) = 



f(L/£(W)), one can replot the data of Fig. ^ by rescaling L. We numerically chose values 
of £ by shifting the logarithm of the length separately for each W. An overlap between data 
for adjacent values of W allows to fit most of the points onto a common curve with two 
branches, one growing (W > W c ) and one decaying (W < W c ) (Fig. [|). Thus we find a 
function £,(W) up to an arbitrary factor £ G . This arbitrariness can be avoided. The points 
for the maximum disorder W = 20 in Fig. ||] except of the sizes L = 6 and 8 correspond to 
the localized regime £ ^ L. Supposing that the statistics of the levels of localized states is 
not sensitive to the dimensionality of system, one can match £ for the above disorder to the 
relevant values of the localization length calculated by the TM-method in the quasi-lD case 
| T5fl . The function £(W), as found by this procedure, is plotted in the inset of Fig. |3|. Our 



values of £ coincide with those found by the TM-method fI7||T6|j . The agreement is slightly 
better for the insulating than for the metallic side. 

The next step is to find the critical exponent. Using the singularity of the correlation 
length near W c 

W -W, 



Z(w) = i 



(3) 



where the scale I is of order of the lattice constant, one can expand the scaling function 
/(L/£) into a power series a(W, L) = a c (L) + A(W — W^L 1 ^, By applying a conventional 
X 2 -criterion for the verification of the linearized scaling hypothesis with several parameters 
to fit the data in Fig. |2| we found v = 1.45 ± 0.08. It is well consistent with the critical 
exponent v ~ 1.53 obtained recently by the TM-method |I8|| . A slightly smaller value of 
v ~ 1.34 for the large-s region of P(s) has been computed in |14| by using the level-statistics 



method similar to ours. The dependence £(W) from (^) with v = 1.45 is shown in Fig. |3| 
for comparison. We note that the critical exponent could be also estimated by fitting the 
data for ^(W) to the power-law function without linearizing of oi(W) near W c . However 
the relative numerical errors become markedly larger similar as in the earlier approach JTS 



The accuracy of v can then be improved by considering a narrower interval of W near the 
critical point or by using larger systems. 



In conclusion, we have studied the integrated level spacing distribution as a scaling vari- 
able in order to detect the critical behavior of the statistics of the one-electron levels at the 
disorder-induced metal-insulator transition. The Anderson model was used to calculate the 
spectra of disordered cubic lattices of sizes varying from 6 3 up to 28 3 sites. Diagonalization 
of the Hamiltonian was performed numerically by means of the Lanczos algorithm espe- 
cially optimized for sparse symmetric matrices with a hierarchic structure. We calculated 
the critical distribution P c (s) at the transition and confirmed its size-independence conjec- 



tured in [pLl 1 . Near the critical point the finite-size scaling properties of the probability I(s) 
of neighboring spacings were examined. For the small-s and the large-s parts of P(s) this 
probability exhibits the transition between the Wigner and the Poisson distributions which 
corresponds to the delocalization-localization transition. Using small spacings the critical 
disorder was found, W c = 16.35, for the 'box distribution' of the site energies. Finally, we de- 
termined the disorder dependence of the correlation length £,(W) which is in good agreement 
with previous results obtained by the transfer- matrix method The present calculation 
confirms one-parameter scaling by using a completely independent method, and that the 
critical exponent of the purely disorder-induced metal-insulator transition is v = 1.45 ±0.08. 
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FIGURES 

FIG. 1. Level spacing distributions P(s) for cubes of sizes L 3 = 6 3 and 28 3 for various values 
of the disorder W near the transition. Continuous curves are the Wigner surmise, Eq. (U), and the 
Poisson distribution, Eq. (Q). 

FIG. 2. Scaling variables a(W) and ~/(W) near W c for different L = 6(0), 8(+), 12(D), 16(x), 
20(A), 24(*), 28(o). 

FIG. 3. Scaling variable a as a function of L/^(W), showing one-parameter scaling. 
W = 12(0), 14(+), 15(D), 16(x), 16.5(A), 17(*), 18(o), and 20(«). Inset: correlation length 
£ as function of the disorder W. Present method (•), TM-method (O) |i~5| ; (+) p6| . Continuous 
curve: plot of Eq. @ with W c = 16.35, v = 1.45, and I = 1. 
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